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^ ! Development of a BEM Solver using Exact Expressions for Computing 

^^' the Influence of Singularity Distributions 
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^ , Summary 

I . Closed form expressions for three-dimensional potential and force field due to 

'^' singularities distributed on a finite flat surface have been used to develop a fast and 

accurate BEM solver. The exact expressions have been investigated in detail and 
found to be valid throughout the complete physical domain. Thus, it has been pos- 
^ ] sible to precisely simulate the complicated characteristics of the potential and force 

O ■ field in the near-field domain. The BEM solver has been used to compute the ca- 

^ ■ pacitance of a unit square plate and a unit cube. The obtained values have compared 

(— I i well with very accurate values of the capacitances available in the literature. 

Introduction 

. The boundary element method is based on the numerical implementation of 

^ , boundary integral equations based on the Green's formula. In order to carry out the 

cn I implementation, the boundaries are generally segmented, and these boundary ele- 

2 ' ments are endowed with distribution of singularities such as source, sink, dipoles 

(^ . and vortices. The singularity strengths are obtained through the satisfaction of 

^ I boundary conditions (Dirichlet, Neumann or Robin) which calls for the computa- 

tion of the influence of the singularities at the points where the boundary conditions 
are being satisfied. Once the singularity strengths are known, physical properties 
at any point within the physical domain can be easily estimated. Thus, accurate 
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distributed on a surface is of crucial importance. 



In this work, we have presented a BEM solver based on closed form expres- 
k>( , sions of potential and force field due to a uniform distribution of source on a flat 

V^ ' surface. In order to validate the expressions, the potential and force field in the 

near-field domain have been thoroughly investigated. In the process, the sharp 
changes and discontinuities which characterize the near-field domain have been 
easily reproduced. Since the expressions are analytic and valid for the complete 
physical domain, and no approximations regarding the size or shape of the singular 
surface have been made during their derivation, their application is not limited by 
the proximity of other singular surfaces or their curvature. Moreover, since both 
potential and force fields are evaluated using exact expressions, boundary condi- 
tions of any type, namely, Dirichlet (potential), Neumann (gradient of potential) or 
Robin (mixed) can be seamlessly handled. 
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As an application of the new solver, we have computed the capacitances of a 
unit square plate and a unit cube to very high precision. These problems have been 
considered to be two major unsolved problems of the electrostatic theory im-l|5|. 
The capacitance values estimated by the present method have been compared with 
very accurate results available in the literature (using BEM and other methods). 
The comparison testifies to the accuracy of the approach and the developed solver. 

Theory 

The expression for potential (V) at a point {X,Y, Z) in free space due to uniform 
source distributed on a rectangular flat surface having corners situated at {x\,z\) and 
{x2,Z2) can be shown to be a multiple of 

V(X,Y,Z)= / , (1) 

where the value of the multiple depends upon the strength of the source and other 
physical considerations. Here, it has been assumed that the origin of the coordinate 
system lies on the surface plane (XZ). The closed form expression for V{X,Y,Z) is 
as follows: 
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where 



Ri=Y^ + {Z-Zif\R2 = Y^ + {Z-Z2? 
h = {X-xi) \Y\;l2 = {X-X2) \Y\\Si =sign{z\ -Z);S2 = sign{z2-Z) 

Similarly, the force components are multiples of 

F.(X,y,Z) = in f^^-j^-^'ll - in f^^^^^l (3) 
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In Eq.(4), C is a constant of integration as follows: 

{0 if outside the extent of the flat surface 

2 71 if inside the extent of the surface and Y > 

—271 if inside the extent of the surface and Y < 

A BEM solver has been developed based on Eqs. (2)-(5) which has been used to 
compute the capacitances of a unit square plate and a unit cube. It is well known 
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Figure 1 : Comparison of potential distribution along X axis 



that the charge densities near the edges of these bodies are much higher than those 
far away from the edges. In order to minimize the errors arising out of the assump- 
tion of uniform charge density on each panel, we have progressively reduced the 
segment size close to the edges using a simple polynomial expression as used in 
El- 
Results 

In order to establish the accuracy of the exact expressions, we have compared 
the potential and electric field distributions computed using the new expressions 
with those computed by assuming a varying degree of discretization of the given 
surface, where each of the discrete element is assumed to have its charge concen- 
trated at its centroid. The flat surface has been assumed to be a square (1cm x 1cm) 
and length scale up to lO/^m has been resolved. The strength of the source on 
the flat surface has been assumed to be unity. In FigurefQ we have presented a 
comparison among results obtained by the exact expressions and those obtained by 
discretizing the flat surface having a single element, having 10x10 elements and 
having 100 x 100 elements. The maximum discretization 100 x 100 discretization 
apparently yields good results even close to the origin. However, oscillations in the 
potential value are visible with sufficient close-up. Similar remarks are true for Ey, 
but for Ex (Figur^SJi, even the highest discretization leads to significant amount of 
oscillation in the estimate. From these figures, we can also conclude that the exact 
expressions reproduce the correct features of the fields even in the most difficult 
situations. The figures vividly represent the error incurred in modeling distributed 
sources using the point analogy which is one of the most serious approximations of 
the BEM. 

In Table ^ we have presented a comparison of the values of capacitances for 
a unit square plate and a unit cube as calculated by |[Q-||5l and our estimations. 
Although it is difficult to comment regarding which is the best result among the 
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Figure 2: Comparison of E^ distribution along X axis 



published ones, it is clear from the table that the present solver indeed lead to very 
accurate results which are well within the acceptable range. Finally, in figure |3l the 
variation of charge density on the top surface of the unit cube has been presented. 
The sharp increase in charge density near the edges and corners is quite apparent 
from the figure. 

Conclusions 

A fast and accurate BEM solver has been developed using exact expressions 
for potential and force field due to a uniform source distribution on a fiat surface. 
The expressions have been validated and found to yield very accurate results in the 
complete physical domain. Of special importance is their ability to reproduce the 
complicated field structure in the near-field region. The errors incurred in assuming 
discrete point sources to represent a continuous distribution have been illustrated. 
Accurate estimates of the capacitance of a unit square plate and that of a unit cube 
have been made using the new solver. Comparison of the obtained results with 
very accurate results available in the literature has confirmed the accuracy of the 
approach and of the solver. 
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Method 


Plate 


Cube 


m 


Surface Charge 


0.3607 


- 


m 


Refined Surface Charge 
and Extrapolation 


0.3667892± 1.1x10"" 


0.6606747 ±5 X 10"' 


m 


Refined Boundary Element 
and Extrapolation 


0.3667874±lxl0-' 


0.6606785±6xlO-' 


m 


Numerical Path Integration 


0.36684 


0.66069 


m 


Random Walk on the Boundary 


- 


0.6606780±2.7xl0-' 


This work 


Boundary Element with 
Exact Expression for Potential 


0.3667869 


0.6606746 



Table 1 : Comparison of capacitance values 
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Figure 3: Variation of charge density on the top surface of the cube 



